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ABSTRACT 

Usual supernova remnants have either ionizing plasma or plasma in coUisional ionization equilib- 
rium, i.e., the ionization temperature is lower than or equal to the electron temperature. However, 
the existence of recombining supernova remnants, i.e., supernova remnants with the ionization tem- 
perature higher than the electron temperature, is recently confirmed. One suggested way to have 
recombining plasma in a supernova remnant is to have a dense circumstellar medium at the time of 
the supernova explosion. If the circumstellar medium is dense enough, coUisional ionization equilib- 
rium can be established in the early stage of the evolution of the supernova remnant and subsequent 
adiabatic cooling which occurs after the shock wave gets out of the dense circumstellar medium makes 
the electron temperature lower than the ionization temperature. We study the circumstellar medium 
around several supernova progenitors and show which supernova progenitors can have a circumstellar 
medium which is dense enough to establish coUisional ionization equilibrium soon after the explo- 
sion. We find that the circumstellar medium around red supergiants (especially massive ones) and the 
circumstellar medium which is dense enough to make Type Iln supernovae can establish coUisional 
ionization equilibrium soon after the explosion and can evolve to recombining supernova remnants. 
Wolf-Rayet stars and white dwarfs have possibility to be recombining supernova remnants but the 
fraction is expected to be very small. As the occurrence rate of the explosions of red supergiants 
is much higher than that of Type Iln supernovae, the major progenitors of recombining supernova 
remnants are likely to be red supergiants. 

Subject headings: ISM: supernova remnants — supernovae: general 



1. INTRODUCTION 

X-ray observations by the Suzaku satellite are con- 
firming t he existence of rec ombining s upernova rernnants 
(SNRs) (lYamaguchi et al][20 09. 201 2i:lOzawa et al.|[200l 
lOhnishi et al.l 120111: iSawada fc Kovamal 120121) ~Recom- 
bining SNRs are SNRs in which the ionization temper- 
ature is higher than the electron temperature. The for- 
ward shock wave emerged at the time of a supernova 
(SN) explosion propagates in the interstellar medium 
(ISM). As the typical density of the ISM is very small 
(rie ~ 1 cm~^ or less where Ue is the electron number den- 
sity) , the timescale to reach coUisional ionization equilib- 
rium (CIE) in t he shocked I SM is typically ^ 10'^ years 
or longer (e.g., iMasail [r984[ ) . Electrons heated by the 
Coulomb interaction with ions in the shocked ISM col- 
lisionally excite ions and reach CIE with this timescale. 
Thus, young SNRs before CIE are supposed to be ion- 
izing SNRs in which the electron temperature is higher 
than the ionization temperature and evolve to SNRs in 
CIE. Most of SNRs are known to be in either the ion- 
izing stage or CIE (e.g., iKawasaki et al.l [20051 ) . In this 
simple picture, SNRs cannot be recombining SNRs and 
the confirmation of the recombining SNRs challenges the 
current understanding of the evolution of SNRs. 

There are several suggested mechanisms to make re- 
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comb ining plasma in SNRs (see, e.g., lYamaguchi et al.1 
120121 and references therein). The existence of a dense 
circumstellar medium (CSM) is one possible way to ex- 
plain the recombin i ng SNRs (e.g., | Itoh fc Masail 119891 : 
IShimizuet aT1[20TTI : IZhou et al.ll201ll) . If a deme CSM 
is around a SN, CIE can be achieved in much shorter 
timescale [~ 10''^/(ne/l cm~^) years]. When the shock 
wave reaches the outer edge of the dense CSM, the 
shocked CSM suddenly expands adiabatically and the 
electron temperature suddenly becomes low and the 
plasma starts to recombine. 

Although the existence of the dense CSM at the time of 
the SN explosion has been suggested as a possible mech- 
anism to realize recombining SNRs, we still do not have 
a clear picture about possible SN progenitors that can 
have a CSM which is dense enough to make recombining 
SNRs. Here, in this Letter, we look into the properties of 
the CSM around SN progenitors at the time of SN explo- 
sions and investigate the progenitors which can evolve to 
recombining SNRs. We focus on massive star progenitors 
because two recombining SNRs, IC 443 and W49B, are 
clearly associated with massive star forming regions (e.g., 
lYamaguchi et al.ll2012D but we also investigate possible 
channels for white dwarfs to be recombining SNRs. 

2. POSSIBLE PROGENITORS 

2.1. Red Supergiants & Wolf-Rayet Stars 

Red supergiants (RSGs) and Wolf-Rayet (WR) stars 
are progenitors of core-collapse SNe. Because of their 
high luminosities, they lose their mass before their ex- 
plosions. Thus, RSGs and WR stars explode inside the 
CSM created by the preceding stellar evolution. If we 
assume that the CSM is from a steady wind with the 
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Fig. 1. — Evolution of the electron number density in the shocked 
CSM from RSG winds. Vs = 10, 000 km s"! and v^ = 10 km s"! is 
fixed and lines with several mass-loss rates are shown. The distance 
in the horizontal axis is the distance which can be reached by the 
forward shock with Vs = 10, 000 km s~^ in the corresponding time. 
If we assume that the progenitor is in the RSG stage for 10 years 
before the explosion, the CSM can reach 1 pc with Vn; = 10 km s~^ . 
The dashed line follows Uest = lO^'^ cm~^ sec and corresponds to 
the minimum density required to reach CIE at the time. Shocked 
CSM above this line is presumed to be at CIE. 



velocity Vw and the mass-loss rate M, the wind density 
Pw becomes 

M 



Airr^v,, 
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where r is the radius. 

If the star inside the CSM explodes, a forward shock 
propagates in the CSM. Assuming that the adiabatic in- 
dex of the system is 3/5, the density ps of the shocked 
CSM just behind the forward shock becomes ps = 4^^,. 
As the forward shock propagates in the CSM, it is decel- 
erated, especially if the CSM is dense. However, the mass 
of the CSM swept up by the forward shock is still small 
compared to the progenitor mass in the early epochs wc 
arc interested in and we assume that it is freely expand- 
ing with the velocity Vg for simplicity. Note that the 
deceleration makes the time of the interaction between 
the shock wave and the CSM longer and CIE can be 
achieved easier with the deceleration. The typical Vs 
of standard SN explosions is Vg ~ 10,000 km s~^ (e.g. 



Suzuki fc Nomotolll995HFransson et al.lll996l:[DwarkadaJ 
20051I2007D . The location of the forward shock at the 
time t after the explosion is r = Vgt and ps can be ex- 
pressed as 

Although Equation ^ is the evolution of the density just 
behind the shock, the remaining entire shocked CSM has 
similar densities when the shock is trave ling in the den- 
sity structure cl ose to p^ j oc r~^ (see, e.g.,|C hcvalicr"198_ 
Suzuki i^^Nomoto 1995i : iFransson et al.|[l996: Dwarkada; . 
20051 120071 ) and we assume that ps is a typical value in 
the shocked CSM. The actual densities in the shocked 
CSM are slightly higher than ps. 

Since the wind properties of RSGs and WR stars differ, 
we consider two cases separately. 



2.1.1. Red 

The typical mass-loss rate and wind velocity of RSGs 
are ^ 10~^ M(7i yr~^ and^ 10 km s^^, respectively (e.g., 
IMauron fc Josselml 1201 If ) and they are consistent with 
those estimated from the obser vations of SN explosi ons 
from RSGs (Type IIP SNe, e.g.. lChevalier et al.l[2006l) . If 
we assume that the RSG wind has the solar metallicity 
and H and He in the wind are fully ionized when the 
forward shock passes, ps = 2.0 x I0~^^nes where Ues is 
the electron number density in the shocked CSM. From 
Equation ^ , the time evolution of Ues is 
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where A/_ 5 is Af scaled by 10 ^ A/q yr ^, ■i;s,9 is w^ scaled 
by 10,000 km s~^, and u^.g is v^ scaled by 10 km s~^. 

Electrons and ions in plasma can reach CIE with the 
timescale of 

Uest ^ 10^^ cm"^ sec, (4) 

(e.g.. lMasall984l:ISmith fc Hughesl20Tol) . Note that only 
ions are presumed to be heated by the forward shock 
and electrons are heated up by the subsequent Coulomb 
interaction between ions an d electrons. T he timescale of 
the electron heating is (e.g.. lMasailll994[ ) 
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where p is the mean molecular weight and In A is the 
Coulomb logarithm. Although the timescale of temper- 
ature equilibrium is a few orders of magnitudes longer 
than that of CIE, the electron temperature can reach 
about 10 % of the ion t emperature {^ 10^ K) in the CIE 
timescale ()Masailll994[ ) and becomes high enough to ex- 
plain the ionization temperature of recombining SNRs. 
Recombining plasma in SNRs can appear if electrons c ool 
down after CIE is achieved (e.g.. lItoh fc Masail [19891 ) . 

Figure [1] shows the comparison of the evolution of the 
typical density in the shocked CSM (Equation (U))) and 
the CIE timescale (Equation ^). CIE can be achieved 
at early epochs of SNRs from RSGs with the typical 
mass-loss rate if we take into account the existence of the 
CSM. This is contrary to the general belief that it takes 
much time to be CIE because SNRs evolve in ISM. Note 
that Vs in the early time is presumed to b e higher than 
the value assumed in Figure [1] (e.g., iDw arkadas 200^) 
and this effect can make the evolution of the electron 
number density faster. In addition, the mass of the re- 
combining plasma estimated from Figure [T] in the case 
of the standard mass loss is ~ 3 x 10~^ Mq and rather 
small. Massive RSGs, yellow supergiants, or RSGs in 
binary systems can have higher mass-loss rates than less 
mas sive RSGs esp ecially just before their explosions (see, 
e.g.. lGeorgvll2012l and references therein) and they are 
more likely to become recombining SNRs among RSGs. 

The early X-ray observations of Type lib SN 1993J 
whose prog enitor is an RSG in a binary system (e.g., 
iMaund et a l. 2004) revealed the existence o f the CIE 
plasm a at a few days since the explosion (e.g.. lUno et al.l 
l2002f ) and the progenitor's mas s-loss rate is suggested to 
be - 5 x 10-5 Mr;, yr-i (e.g.. ISuzuki fc Nomotol [l995l: 
IFransson et al.|[i996f ). This is consistent with our esti- 
mate and explosions of RSGs can establish CIE in the 
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Fig. 2. — Evolution of the electron number density in the 
shocked CSM from WR star winds. Vs = 10,000 km s~^ and 
Vm = 1000 km s~^ is fixed and lines with several mass-loss rates 
are shown. The distance in the horizontal axis is the distance which 
can be reached by the forward shock with Vs = 10, 000 km s^^ in 
the corresponding time. If we assume that the progenitor is in 
the WR stage for 10^ years before the explosion, the CSM can 
reach 100 pc with d™ = 1000 km s~^. The dashed line follows 
nest = 10 cm~^ sec and corresponds to the minimum density 
required to reach CIE at the time. Shocked CSM above this line 
is presumed to be at CIE. 



early epochs and evolve to recombining SNRs. 

2.1.2. Wolf-Rayet Stars 

The typical mass-loss rate and wind velocity of WR 
stars are ~ 10~^ Mp, yi~^ and ~ 1000 km s~^, respec- 
tively (e.g.. lCrowther|[2007l ). If we assume that the wind 
from WR stars is composed of 50 % of carbon and 50 
% of oxygen and they are fully ionized when the CSM 
is shocked, Ps ~ 3.3 x 10~^*nes- With these values, the 
following equation is obtained from Equation ([2]): 
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where v^fi is Vw scaled by 1000 km s^^. 

The evolution of rig^ is compared to the CIE timcscale 
TT-esi ^ 10^^ cm~'^ sec in Figure [5J With the canoni- 
cal mass-loss rate 10~^ Mq yr~^, the evolution of the 
shocked CSM can be comparable to the CIE timescale 
just after the explosion. However, Vs is presumed to be 
larger than Vs = 10,000 km ^~^ in the early epochs and 
the evolution of ngs can be faster, as noted in the pre- 
vious section. Contrary to the case of RSGs, CIE can 
be only established at the very early epochs with the op- 
timistic Vs and it is likely that CIE is not achieved at 
early phases of the typical explosions of WR stars with 
the canonical mass-loss history. This is because the typ- 
ical wind velocity is about 100 times larger than the typ- 
ical RSG wind velocity and the CSM becomes thin much 
faster. Although recent radio observations of explosions 
of WR stars (Type Ibc SNe) are revealing the existence 
of W R stars with high ma ss-loss rates (~ 10~* Mq yr~^, 
e.g., IWellons et al .1120 121) . the amount of recombining 
plasma is very small (^ 10""' Mq, Figure [2]) even if such 
a high mass-loss rate is maintained for the entire WR 
phase. Thus, WR stars may be difficult to have dense 
CSM which is massive enough to make them rccombining 
SNRs. Since some elements can reach the CIE in smaller 



timescalcs than rif^s^ ~ 10^^ cm ^ sec at the typical tem - 
pcrature in the shocked CSM (jSmith fc Hughes! I2010l) . 
at least some elements may reach CIE. In addition, it 
is also known that some WR stars experience explosive 
mass loss just before their explosions which can eject 
massive CSM, as is indicated by the proge nitor of Type 
lb SN 2006JC (e.g., iPastorello et~all l2007l) . Explosions 
of this kind of WR stars can also result in recombining 
SNRs but they are also expected to be rare. 

2.2. Type Iln Supernova Progenitors 

Type Iln SNe are SNe which show narrow spectral 
lines, especia lly of hydrogen, in optical spectra (e.g., 
lSchlegellll990[) . Their spectral features can be explained 
by the existence of the dense CSM wit h the Thom- 
son optical depth t t larger than 1 (e.g., iChugail [20011 : 
iDessart et alJ 120091 ). They can become luminous in X- 
ray and radio and some of them become luminous even 
in optical. The high luminosities of Type Iln SNe can 
be naturally explained by the i nteraction between SN 
ejecta and its den se CSM (e.g.. iDwarkadas et al.ll20l"ol : 
iChugai et"ani2004[ ) and the progenitors of Type Iln SNe 
are a plausible candidate of the origin of recombining 
SNRs. 

We assume that Type Iln SNe have a dense CSM with 
Tt = arriew^R '^ 1, where ctt is the Thomson cross sec- 
tion, n^u, is the mean CSM electron number density, and 
AR is the CSM length. We use the mean density riew 
because the mass loss of Type Iln SN progenitors just be- 
fore their explosions are revealed to be no n-steady from 
X-ray observations (jDwarkadas fc Gruszko 2012). The 
shock wave with the velocity Vs can propagate through 
the CSM with ts = AR/vs. As CSM or ISM with 
much lower density exists outside the dense CSM and 
the shocked CSM is rarefied after the shock wave goes 
out of the dense CSM, recombining plasma can be eas- 
ily synthesized once the CIE is achieved in the shocked 
CSM. Assuming Ties = 4neMi and the solar metallicity, 
the typical timescale before the rarefaction is 
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This is much larger than the timescale required to achieve 
the CIE, riet ~ 10^^ cm~'^ sec. In reality, Vs can be 
smaller than 10,000 km s^^ because of the deceleration 
by the dense CSM but this makes ts longer. One caveat 
is that we use a constant mean electron density ne™ to 
estimate the density evolution. If the density declines 
very steeply, this assumption can be very crude and Type 
Iln SNe from very steep CSM may not end up with re- 
combining SNRs. Nonethele ss, many Type Iln S Ne have 
flat density CSM fJDwarka das fc Gruszkoll2012n and we 
presume that most of Type Iln SNe can end up with re- 
combining SNRs. However, as explosions of RSGs (Type 
II SN e) occurs muc h more frequently than Type Iln SNe 
fe.g.. lLi et al.|[20Tl|) . the major progenitors of recombin- 
ing SNRs are likely to be RSGs. 

Unfortunately, the progenitors of Type I ln SN e are 
not understood well. I Gal- Yam fc Leonard! ()2009! ) have 
confirmed that the progenitor of a Type Iln SN, 
SN 2005gl, is a lumiiious b lue variable (LBV, e.g., 
!Humphrevs fc DavidsonI !1994| ). As LBVs originate 
from very massive stars (more than Mzams '^ 40Mq 
where Mzams is the zero-age main-sequence mass, e.g.. 
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iCrowtherl ()2007[ )) and the observational rate of Type Iln 
SNe is also consist ent with the mass range of LBVs (e.g., 
iSmith et al.ll201lf ). Type Iln SNe are suggested to come 
mainly from these very massive stars. However, it is the- 
oretically considered that LBVs arc in an evolutionary 
stage in which very massive stars evolve to WR stars 
and LBVs do not explode. It is also possible that the 
fast wind from a WR star collides to the slowly mov- 
ing wind from its previous RSG stage and a dense shell 
which is enough to be a Type Iln SN is created by 
the interaction (e.g.. iDwarkadas et al.l 120101 ). Another 
possible Type Iln SN progenitor is a super-asymptotic 
giant branch (AGB) star with Mzams ~ 8 M©. An 
0-|-Ne-|-Mg core at the cent er of the super- AGB wind can 
be an electron-capture SN ()Nomotdll984[ ). The progen- 
itor of Type Iln SN 2008S is found to be around 10 Mq 
(IPrieto et a l. 20081 and may belong to this class (e.g., 
iBotticella ct al. 2009l) but there also exists an argumen t 
that SN 2008S may not be an SN (|Smith et all 120091) . 
Finally, we note that some o f Type Ic superlumi nous 
SNe recently discovered (e.g., iQuimbv et al.l [20Tl|) are 
related to the inter action of a dense C+ O-rich CSM and 
SN el ecta (Blinnik ov fc Sorok ina 2010; Moriva fc Maedal 
|2012| ) and they can also be a progenitor of recombining 
SNRs. However, Type Ic superlum inous SNe preferen- 
tially appear in m etal-poor galaxies (jQuimbv et al.|[20Tll : 
iNeill et all 120 111 ) and the occurrence rate is also quite 
small. Thus, recombining SNRs currently observed in 
our Galaxy seem irrelevant to them. 

2.3. White Dwarfs 

Although most of the recombining SNRs currently dis- 
cover ed are likely to origina te from core-collapse SNe 
(e.g., lYamaguchi et al.l I2012D . Type la SNe can also 
evolve to recombining SNRs although it is expected to 
be quite rare. Type la SNe are explosions of white 
dwarfs. There are two major suggested paths for white 
dwarfs to exp lode: single degenerate (SD) channel (e.g. 
Nomotolll982D and dou ble degenerate (DD) channel (e.g. 



Iben fc Tutukovlll984l ). In the SD scenario, a white dwarf 



is in a binary system with a main-sequence star and the 
mass of the companion accretes to the white dwarf. The 
white dwarf explodes when its mass gets close to the 
Chandrasekhar mass limit. On the other hand, the DD 
scenario suggests that Type la SNe are caused by the 
merger of two white dwarfs in a binary system. The 
main channel of Type la SNe is still unknown. 

In the SD scenario, the exploding white dwarf is sur- 
rounded by the accr eting materials with a typical rate 
of ~ 10-"^ Mq yr-i (|Nomotolll982l ) but the rate is pre- 
sumed to be too small to make the recombining SNR. 
Mass loss from the companion star can also make CSM 
around the progenitor but the companion is likely to be a 
less evolved red giant with too small mass-l oss rates to es- 
tablis h CIE (~ 10~^ Mq yr~^ or less, e.g.. lHachisu et all 
Il999| ). However, there are rare ways to make the mass- 
loss rate of the syste m high during the binary evolution 
(jHachisu et al.l 120081 ) and some Type la SNe are actu- 
ally suggested to be a hybrid of Type la and Type Iln, 



i.e.. Type la SNe exploded in a CSM as dense as those 
discu ssed in Section 12.2! (e.g., SN 2002ic, iHamuv et al.l 
I2003D . Thus, it is possible that a Type la SN from the 
SD scenario evolve to a recombining SNR but the number 
is expected to be very small. 

In the DD scenario, we do not expect the CSM from 
the progenitor system because two binary stars are white 
dwarfs. However, stripped materials at the time of the 
merger are sugg ested to remain w hen the merged white 
dwarf explodes (jFrver et al.ll2010l ). These materials are 
quit e dense (uew > 10\^ cm~^ within r = Rq, see Figure 
5 of iFrver et all (|2010[) 1 and Type la SNe exploded in 
such environment can reach CIE and may end up with 
recombining SNRs. However, such a den se envelope is 
not ob tained in a similar DD simulation of lPakmor et al.l 
()2012D . The fact that wc do not sec recombining SNRs 
of Type l a SNe may already suggest that the model ob- 
tained by iFrver et all ()2010r) is not the major path to be 
a Type la SN. Because of the uncertainty in the theo- 
retical prediction of Type la SNe from the DD scenario, 
we still cannot exclude the possibility that Type la SNe 
from the DD scenario can be recombining SNRs. 

To sum up, although all the recombining SNRs cur- 
rently discovered are likely from core-collapse SNe, Type 
la SNe from both the SD and DD channels have possi- 
bility to become recombining SNRs. When the detailed 
theoretical predictions are fixed, like the existence of the 
dense envelope in the DD channel, recombining SNRs 
may be able to be a probe to indicate the progenitor 
system of Type la SNe. 

3. CONCLUSIONS 

We have investigated the possible progenitors of recom- 
bining SNRs. If a CSM which is dense enough to estab- 
lish CIE in the early epochs of the SNR evolution exists 
around a progenitor, the plasma in the shocked CSM can 
be overionized and the SNR can become a recombining 
SNR. RSGs, especially massive ones, and Type Iln SN 
progenitors can have the CSM which is dense enough to 
establish CIE at the early stage of their explosions and 
can evolve to recombining SNRs. As explosions of RSGs 
(Type II SNe) occurs much more frequently than Type 
Iln SNe, the major progenitors of recombining SNRs are 
likely to be RSGs. 

WR stars and white dwarfs are difficult to make re- 
combining SNRs with their standard mass-loss histories 
but they are suggested to have mechanisms to enhance 
their mass-loss rates and they can be recombining SNRs 
if such mechanisms enhance their mass-loss rates. How- 
ever, these mechanisms arc presumed to work on a small 
fraction of these stars and such progenitors are expected 
to be a minor way to have recombining SNRs. 
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